Ferromagnetism in the Hubbard model with orbital degeneracy in infinite dimensions 
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We study the ferromagnetism due to orbital degeneracy in the Hubbard model in infinite dimensions. 
The model contains the intraorbital repulsion U , the interorbital repulsion U' , the exchange J (Hund 
coupling) and the pair hopping J', where all of them originate from the on-site Coulomb interaction. 
The ground state of the effective one-site problem was obtained by exact diagonalizations. At the 
1/4-filling, we found two insulating phases; one is a ferromagnetic phase with alternating orbital 
order and the other is antiferromagnetic with uniform orbital order. If electrons are doped into the 
1/4-filling, the ferromagnetic phase still survives and becomes metallic, while the antiferromagnetic 
phase disappears. This result indicates that the double-exchange mechanism is relevant to stabilize 
metallic ferromagnetism in infinite dimensions. 



PACS numbers: 75.10.Lp. 

Many investigations have been done on the Hubbard 
model, to clarify whether the Coulomb interaction can re- 
alize ferromagnetism in itinerant electron systems. They 
revealed that the simple single-band Hubbard model on 
the hypercubic lattice does not easily show ferromag- 
netism and that some additional properties are necessary 
to stabilize ferromagnetism. There are several propos- 
als for such properties. One is the special lattice struc- 
ture (or band structure) which favors ferromagnetism. 
For example, ferromagnetism on the Jop lattices and 
of the flat (or nearly flat)-band modeloa are discussed 
in this context. n Another is the existence of the or- 
bital degeneracyoH It was argued that the on-site Hund 
coupling between electrons in degenerate orbitals cre- 
ates indirect ferromagnetic coupling between conduction 
electrons. □ 

Numerical study of strongly correlated electron sys- 
tems is quite difficult in higher dimensions than one. 
In the infinite-dimensional (d = oo) limit, however, 
it is manageable with the aid of recently developed 
technique!] In this limit we can treat quantum fluctua- 
tions completely, by taking local interactions into account 
exactly, and we can neglect spatial fluctuations. Monte 
Carlo simulations of the d = oo single-band Hubbard 
model revealed that any ferromagnetic phase does not 
appear on the hypercubic latticed and, on the other hand, 
metallic ferromagnetism occurs on the fcc-type lattice.H 
These two results clearly show the lattice-structure de- 
pendence for the appearance of the ferromagnetism. The 
purpose of this paper is to clarify whether the itinerant 
ferromagnetism appears due to orbital degeneracy in in- 
finite dimensions. 

We study the Hubbard model with doubly degener- 
ate orbitals at each site on the hypercubic lattice. The 
Hamiltonian is given by 
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where Ci ma (c| mo .) denotes the annihilation (creation) op- 
erator of the electron at site i with orbital m(=l or 2) 
and spin a. The number operator is denoted by ft,- mCT . 
All the interaction terms in Eq. ([!]) originate from the 
Coulomb interaction between electrons at the same site. 
We treat every interaction exactly. Generally interaction 
parameters satisfy the relation U > U' > J ~ J'. 
symmetry argument leads to the relation U = U' + 2 J 
Hopping terms are scaled with the dimensionality d so 
that the one-electron density of states (DOS) behaves as 
D(e) = exp(— e 2 /t 2 )/ty/ir in the limit d — > oo. 

The double-exchange mechanism due to strong Hund 
coupling J was proposed to produce itinerant ferromag- 
netisirL.for n > 1, where n denotes electron density per 
site.tatl At the 1/4-filling (n — 1) the ground state is ex- 
pected to be insulating in the strong-coupling limit. In 
this case the effective interaction due to virtual hoppings 
of electrons favor a ferrjOmacnctic ground state with al- 
ternating orbital order.EirE-3 In one dimension, the 
romagnetism was found numerically near n = 1,lJ 
and shown rigorously for general n in the strong-coupling 
limit .EJ'EZl For higher dimensions, it is not yet clear 
whether ferromagnetic order appears. Inagaki and Kubo 
obtained ferromagnetic phases on the simple cubieJattice 
for J' = using the Hartree-Fock approximation.Eil How- 
ever, it is not quite clear how the predicted metallic ferro- 
magnetism in three dimensions is stable, since the the lo- 
cal quantum fluctuations and the pair-hopping (J') term, 
both of which suppress the magnetic order, are neglected 
in their study. Effects of orbital degeneracy were also 
studied with a d = oo two-band model in the .anisotropic 
(Ising) limit of Hund coupling and J' = 0.E3 O 
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In this paper we study the model (|]) in infinite dimen- 
sions near the quarter filling. We present the one-site ef- 
fective action of the model (Q) and study the ground-state 
properties using the exact diagonalization method. We 
found a metallic ferromagnetic phase for electron doped 
cases, e.g. n — 1.2, which indicates that the double- 
exchange mechanism can stabilize ferromagnetism in in- 
finite dimensions. At the 1/4- filling, the ground state 
shows a metal-insulator transition at a finite parameter 
U . A phase diagram with both ferromagnetic and anti- 
ferromagnetic phases is obtained for n = 1. The phase 
diagram is consistent with a mean-field theory based on 
the perturbational treatment from the strong-coupling 
limit. 

To derive the one-site effective action of the present 
model, we employ the cavity method.EI Taking the trace 
for all degrees of freedom except for those on site 0, we 
obtain the effective action, 



exp[ 



SeS (Cqtoct' c 0ma)] 



~Z~ 
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where the partition function of the model is written as a 
functional integral over Grassman variables. Generally, 
S e ff is an infinite series of multiples of Grassman variables 
at various times. In the d = oo limit, the action becomes 
quite simple as 



S, 



cff 



EE- 1 

m=l,2 u 



dTH- mt (cl ma {T),CQ ma (T)), 



(3) 



where H- mt denotes the local interaction on the site 0, 
and Q{) ma satisfies the relation 



(4) 



To derive Eqs. (|) and (@), we have used the facts that 
the self-energy is independent of momentum and that 
(c\cj) ~ (I / \fd) <[% ~ i\ , which applies for large d. 

Here we briefly explain the method for calculating the 
ground state. This one-site effective action can be ap- 
proximated by the following two-channel impurity model: 



m=l,2 tp=2,n s 

p=2,n s ) 
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where a pmcr (c ma ) denotes the annihilation operator of 
the conduction (impurity) electron with channel m and 
spin a. The effective action of the impurity site has the 
same form as Eq. (||) with 
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instead of Qoma-(iw n ). Conduction electrons have n s — 1 
orbitals. The effective action (||) can be well approxi- 
mated with the impurity model (^), if we increase n s 
and select the most suitable parameters {ep m a,Vpma) so 
that (?™p satisfies the self-consistency condition (Eh. We 
numerically diagonalize the impurity model (JsJ) with the 
Lanczos method and obtain the ground state under the 
self-consistency condition (|J). It has been shown that 
this method works well and the ground state can be ac- 
curately described even with small number n s = 5 or 



6.0 In spatially uniform ground states the Green func- 
tion is related to the self-energy with G mcr (k, zcj n ) _1 = 
(Go) -1 — £ — £fe, where Go = l/(iu) n + //) and Ek = 

d~ x / 2 53u=i cos ^m- We also consider states with two- 
sublattice structures. The Green function from I sub- 
lattice to V sublattice, G^ CT (k, iu) n ), has the form 



(Go)- 1 -^ 



(Go)^- E B 



(7) 



We have two effective actions for two sublattices. The 
self-consistency condition is unchanged as QQ ma .(iuj n ) = 
[G^ (T (iw n ) _1 +i;J n (icJ„)] _1 , where g l 0ma (iuj n ) denotes the 
dynamical mean field applied to the sublattice I. 

We study the ground states at the fillings n = 1, 1.2, 
and 0.8, controlling the chemical potential. Parameters 
are set as U = U' + 2 J and J — J' . Both spatially 
uniform states and those with the two-sublattice struc- 
ture are considered. Numerical calculations were done 
for n s = 5 or n s = 6. We mostly studied the system 
with n s — 5 and confirmed phase boundaries using the 
system with n s — 6. We found that n s dependence of 
the ground state is small and the phase diagrams are the 
same between the cases n s = 5 and 6. 

For the 1/4-filling, we found three ground-state phases, 
i.e., paramagnetic, ferromagnetic, and antiferromagnetic 
ones. We obtained the phase diagram shown in Fig. [|. 
Though J will be positive in real systems, we consider 
both positive and negative J to clarify the effects of Hund 
coupling. Near the phase boundary, two (or three) solu- 
tions coexist and their energies cross over. We selected 
the ground state comparing the energies and determined 
the phase diagram. 

For a wide parameter region, we found a ferromag- 
netic phase. This state has an orbital super-lattice struc- 
ture. The pseudospin rf = X] C r( nil < 7 — n i2o-)/2, which 
represents the orbital degree of freedom, is antiferro- 
magnetically ordered. In this phase the compressibility, 
dn/dfi, is always vanishing and hence the ground state 
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FIG. 1. Magnetic phase diagram at the 1/4-filling (n = 1). 
Parameters are fixed as U — U' + 2 J and J = J'. 



is insulating (See Fig. |2|). Thus magnetic transition and 
metal-insulator transition occur at the same coupling pa- 
rameter. The magnetization is hardly reduced from the 
full polarization. For example, the magnetization per 
site(= N e ) is (M)/N e = 0.5000 for U = 18, U' = 10 and 
J = 4, where N e denotes the number of electrons. Rigor- 
ously speaking the spins cannot be fully polarized in the 
d — > oo limit because of the Gaussian DOS.Eil Since the 
reduction is too small, it is invisible within accuracy of 
the present calculation. 

Another insulating phase appears for negative J. In 
this phase the ground state has antiferromagnetic long- 
range order and also uniform orbital order. The sublat- 
tice magnetization per site is very close to 1/2. 

The phase boundaries of the above two phases are well 
understood in terms of the effective Hamiltonian in the 
strong-coupling limit. We start from the ground states 
in the atomic limit, (t = 0) and treat the hopping term 
as a perturbation.E3 Then the second-order perturbation 
leads to the following effective Hamiltoaian for spin op- 
erators Si and pseudospin operators Tjli 2 ]: 
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Clearly, spin interactions have SO(3) symmetry and, un- 
der the conditions U = U' + 2 J and J = J', pseu- 
dospin interactions also have SO (2) rotational symmetry 
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FIG. 2. Chemical potential dependence of the electron den- 
sity at the parameters U = 18, U' — 10 and J = J' = 4. Mag- 
netic order of each state is written by FM, P, Fl and AFI. FM 
and P denote ferromagnetic and paramagnetic metal, respec- 
tively. Fl means ferromagnetic insulator with alternating or- 
bital order and AFI denotes antiferromagnetic insulator with 
no orbital order. 



about the j/-axis. The derivation of Eq. is valid for 
—U'/3 < J < V, where intermediate states have higher 
energy than the unperturbed ground states. The per- 
fect ferromagnetic state with perfect antiferromagnetic 
pseudospin order has the energy — t 2 / (2U' — 2 J) per site, 
while that of the Neel state with perfect ferromagnetic 
pseudospin order is —t 2 U/(2U 2 — 2 J 2 ). Thus the ferro- 
magnetic state has a lower energy for < J < U' and 
the antiferromagnetic state is lower for —U'/3 < J < 0. 
We see in Fig. |l| that the phase boundaries are close to 
the three lines, J = -U'/3, J = and J = U' . The 
present result reflects the fact that the mean-field theory 
is exact in infinite-dimensional localized spin models. 

In order to examine the effects of electron doping to the 
ordered states at the 1/4-filling, we studied the case with 
n = 1.2. At this filling we obtained two metallic phases, 
both of which are spatially uniform, i.e. the ferromagnetic 
phase and the paramagnetic one (see Fig. ||). We could 
not find any antiferromagnetic phase at this filling. The 
transition from the paramagnetic phase to the ferromag- 
netic one is of first order. Both ferromagnetic and para- 
magnetic states exist as solutions of the self-consistency 
equations in a region close to the phase boundary and 
their energies cross over at the phase boundary. 

The ferromagnetic state is metallic, since it has fi- 
nite compressibility (Fig. ||), and it has no orbital or- 
dering. The reduction of magnetization is larger than 
the insulating ferromagnet at the 1/4-filling. For exam- 
ple, (M)/N e = 0.341 for U = 18, U 1 = 10, and J — 4 at 
n = 1.2. The area of the ferromagnetic phase is reduced 
in the phase diagram compared to that of the insulat- 
ing ferromagnetic phase at the 1/4-filling. On the other 
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FIG. 3. Magnetic phase diagram at the filling n = 1.2. 
Parameters are fixed as U — U' + 2 J and J = J'. 



hand, the antiferromagnetic phase completely disappears 
at this filling. This result is consistent with the well- 
known fact that the itinerancy of electrons (or holes) is 
hardly compatible with the antiferromagnetic long-range 
order. 

As a hole-doped case, we study the ground state for 
n = 0.8. At this filling we found only metallic ground 
states which are translationally invariant. We could not 
find any magnetically ordered phase for J < U' < 20. 
Disappearance of antiferromagnetism is understood as 
the result of the itinerancy of the holes. Self-consistent 
solutions that show ferromagnetism also appear for large 
V, e.g. V — 10, but have slightly higher energy than 
that of the paramagnetic states. Since the probability of 
double occupancy is low for n < 1, the double-exchange 
mechanism due to Hund coupling is less effective than 
in the case with n > 1 and hence ferromagnetism may 
disappear. The ferromagnetic ground &t&tf however ap- 
pears even for n < 1 in one dimension J 1 H 1S I This discrep- 
ancy may be understood with the argument that hia 
dimensionality destroys the perfect ferromagnetism,! 
though we need more study to conclude the absence of 
the metallic ferromagnetism at less than 1/4-filling for 
strong coupling. 

In this study, we found the metallic ferromagnetism, 
in which spins are partially polarized, for positive J only 
at more-than 1/4-filling (n = 1.2). This result should be 
compared with the single-band, Hubbard model, which 
does not show ferromagnetism.13 This indicates that the 
double-exchange mechanism due to the Hund-coupling 
favors ferromagnetism and realizes metallic ferromag- 
netism for strong but finite coupling in infinite dimen- 
sions. The ferromagnetism appears in a parameter region 
where < J < U'. We did not search for ground states 
with larger sublattice structures in the present study. 
Hence possibility of incommensurate states still remains. 



Further extension of the study in infinite dimensions and 
also a reliable phase diagram of the metallic ferromag- 
netism in one dimension, which is still not available, will 
be useful to gain insight into the role of the orbital de- 
generacy in the metallic ferromagnetism in two and/or 
three dimensions. 
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